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The charged current lepton production induced by neutrinos in Fe nuclei has been studied. 
The calculations have been done for the quasielastic as well as the inelastic reactions assuming A 
dominance and take into account the effect of Pauli blocking, Fermi motion and the renormalization 
of weak transition strengths in the nuclear medium. The quasielastic production cross section for 
lepton production are found to be strongly reduced due to nuclear effects while there is about 10% 
reduction in the inelastic cross sections in the absence of the final state interactions of the pions. 
The numerical results for the momentum and angular distributions of the leptons averaged over the 
various atmospheric neutrino spectra at the Soudan and Gransasso sites have been presented. The 
effect of nuclear model dependence and the atmospheric flux dependence on the relative yield of 
to e has been studied and discussed. 

PACS numbers: 25.30Pt neutrino scattering - 13.15.+g Neutrino interactions - 23.40.Bw Weak-interaction 
and lepton (including neutrino) aspects - 21.60.Jz Hartee-Fock and random-phase approximations 



I. INTRODUCTION 

The study of neutrino physics with atmospheric neu- 
trinos has a long history with first observations of muons 
produced by atmospheric muon neutrinos in deep un- 
derground laboratories of KGF in India and ERPM in 
South AfricaQ. The indications of some deficit in the 
atmospheric neutrino flux was known to exist from the 
early days of these experiments but the evidence was 
no more than suggestive due to low statistics of the 
experimental data and anticipated uncertainties in the 
flux calculations 2j. The clear evidence of a deficit in 
the atmospheric muon neutrino flux was confirmed later 
when data with better statistics were obtained at 1MB 3], 
Kamiokandepl and SoudanJJ experiments. The most 
likely cause of this deficit is believed to be the phenom- 
ena of neutrino oscillations in which the neutrinos pro- 
duced with muon flavor after passing a certain distance 
through the atmosphere, manifest themselves as a differ- 
ent flavor. The implication of this phenomena of neu- 
trino oscillation is that neutrinos possess a nonzero mass 
pointing towards physics beyond the standard model of 
particle physics. The evidence for neutrino oscillations 
and a nonzero mass for the neutrinos has also been ob- 
tained in the observations made with solar and reactor 
( anti) neutrinos |8( • 

It is well known that, in a two flavor oscillation scenario 
involving muon neutrino, the probability for a muon neu- 
trino with energy E v to remain a muon neutrino after 
propagating a distance L before reaching the detector is 
given by@- 

2„„ 2 fl-27Am 2 (eV 2 )L(km)\ .. 

P - 1 - sin 2 20sin 2 } '— y - '- 1 

V E v (GeV) ) y ' 

where Am 2 = m 2 — m\ is the difference of the squared 
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masses of the two flavor mass eigenstates and 9 is the 
mixing angle between two states. The oscillation param- 
eters Am 2 and the mixing angle 9 are determined by 
various observations made in atmospheric neutrino ex- 
periments. These include the flavor ratios of muon and 
electron flavors, angular and -g distributions of muons 
and electrons produced by atmospheric neutrinos. The 
first claims of seeing neutrino oscillation in atmospheric 
neutrinos came from the measurements of ratio of ra- 
tios i?„ defined as ^/f>* ata from the observations of 
fully contained (FC) events |3j-|j], but there are now 
data available from the angular and -g distribution of 
the atmospheric neutrino induced muon and electron 
events from SK@, MACRO [13 and Soudan[ll| exper- 
iments which confirm the phenomena of neutrino oscil- 
lations. These experiments are consistent with a value 
of Am 2 w 3.2 x lO~ 3 eV 2 and sin 2 29 w 1. The analysis 
of these data assuming a three flavor neutrino oscillation 
phenomenology have also been done by many authors 

The major sources of uncertainty in the theoretical pre- 
diction of the charged leptons of muon and electron flavor 
produced by the atmospheric neutrinos come from the 
uncertainties in the calculation of atmospheric neutrino 
fluxes and neutrino nuclear cross sections. The atmo- 
spheric neutrino fluxes at various experimental sites of 
Kamioka, Soudan and Gransasso have been extensively 
discussed in literature by many aut hors [l3 | - [l6 | . The 
neutrino nuclear cross sections have also been calculated 
for various nuclei by many authors using different nuclear 
models [l7|-[24|]. The aim of the present paper is to study 
the neutrino nuclear cross section in iron nuclei which 
are relevant for the atmospheric neutrino experiments 
performed at Soudan[I3, FREJUSEl and NUSEX|2| 
and planned in future with MINOS j^J, MONOLITH|H 
and INO[2!| detectors. The uncertainty in the nuclear 
production cross section of leptons from iron nuclei by 
the atmospheric neutrinos are discussed. For our nu- 
clear model, we also discuss the uncertainty due to use 
of different neutrino fluxes for the sites of Soudan and 
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Gransasso which are relevant to MINOS, MONOLITH 
and INO detectors p^-p^. 

The momentum and angular distribution of muons 
and electrons relevant to fully contained events produced 
by atmospheric neutrinos in iron nuclei are calculated. 
These leptons of muon and electron flavor characterized 
by track and shower events include the leptons produced 
by quasielastic process as well as the inelastic processes 
induced by charged current interactions. The calcula- 
tions are done in a model which takes into account nu- 
clear effects like Pauli Blocking, Fermi motion effects 
and the effect of renormalization of the weak transi- 
tion strengths in nuclear medium in local density ap- 
proximation. The model has been successfully applied 
to describe various electromagnetic and weak processes 
like photon absorption, electron scattering, muon capture 
and low energy neutrino reactions in nuclei |24j.|30j-|32^. 
The model can be easily applied to calculate the zenith 
angle dependence and the -g distribution for stopping 
and thorough going muon production from iron nuclei 
which is currently under progress. 

The plan of the paper is as follows. In section-II we 
describe the neutrino(antineutrino) quasielastic inclusive 
production of leptons (e~, (x~, e + , from iron nuclei 
for various neutrino energies. In section-Ill we describe 
the energy dependence of the inelastic production of lep- 
tons through A-dominance model and highlight the nu- 
clear effects relevant to the energy of fully contained 
events. In section-IV, we use the atmospheric neutrino 
flux at Soudan and Gransasso sites as determined by vari- 
ous authors and discuss the flux averaged momentum and 
angular dependence of leptons corresponding to different 
flux calculations available at these two sites. 



Formalism 



II. QUASIELASTIC PRODUCTION OF 
LEPTONS 



Quasielastic inclusive production of leptons in nuclei 
induced by neutrinos has been studied by many au- 
thors where nuclear effects have been calculated. Most 
of these calculations have been done either for 16 O 
relevant to 1MB and Kamioka experimentsQ-Qor for 
12 C relevant to LSND and KARMEN experiment p|. 
These calculations generally use direct summation 
method (over man y n uclear excited sta tes) Il7l . clo- 
sure approximation jl8j . Fermi gas model [ljj-[2(|, re l- 
ativistic mean field approximatio n^ ll . continuum ran- 
dom phase approximation fCRPA)|22l| and local density 
approximation piill^ . The calculations for 56 Fe nucleus 
have been reported by Bugaev et al. 1171 i n a shell model 
and in Fermi gas model by Gallagher 34| and Berger et 
al.[3^|. In this section we briefly describe the formal- 
ism and results of our calculations done for quasielastic 
inclusive production of leptons for iron nuclei. 



In local density approximation the neutrino nucleus 
cross section a(E u ) for a neutrino of energy E v scattering 
from a nucleus A(Z, N), is given by 



dp 



a A {E v ) = 2 / dr-^n n (p,r)a N (E u ) 



(2) 



where n n (p, r) is the local occupation number of the ini- 
tial nucleon of momentum p (localized at position r in 
the nucleus) and a N {E v ) is cross section for the scatter- 
ing of neutrino of energy E v from a free nucleon given by 
the expression 



a N {E v ) 



d 3 k' to,, to„ M„ M, 



(2tt) 3 E v E e E n E p 

x Y J Y.\ T \ 2 ^ E -- Ei+E -- E ^ ( 3 ) 

where T is the matrix element for the basic process 

n (k)+n(p)^l-(k')+p{p'), l = e^ (4) 
written as 

T= ^§cos# c u(*/)7 M (l-7s)u(fc) (Jn C ° (5) 

(J M ) CC is the charged current(CC) matrix element of the 
hadronic current defined as 



(jn CC = u(p') FY{q 2 )Y + F^{q 2 )ia^ 



2M 



(6) 



+FX(g 2 ) 7 V] u(p) 



q 2 (q = k — k') is the four momentum transfer square. The 
form factors (q 2 ), F% (q 2 ) and F\(q 2 ) are isovector 
electroweak form factors written as 

FY(q 2 ) = F[(q 2 ) - F^q 2 ), F? (q 2 ) = F?(q 2 ) F?(q% 
FX(q 2 ) = F A (q 2 ) 



where 
Fr(q 2 ) 



(I r. 



(1 4M 2 ) 



G P E n (q 2 )-^G^(q 2 ) 

[gW)-g£V)] 



G p M (q 2 ) = (1 + Vp)GW), G n M (Q 2 ) = VnG p E (q 2 ) 



(7) 



G n E {q 2 ) = {^)^Gl(q 2 %,: U = 



1 — X 1 
x A ™ 4M 2 
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Hp = 1.79, fin = -1.91, M v = 0.84Gey, and A„ = 5.6. 
The isovector axial vector form factor F^(Q 2 ) is given by 



F A (Q 2 



Fa(0) 



(1 



-) 



where M A = l.032GeV; F j4 (Q)=-1.261 

In a nuclear process the neutrons and protons are not 
free and their momenta are constrained by the Pauli 
principle which is implemented in this model by requir- 
ing that for neutrino reactions initial nucleon momentum 
p < pf„ and final nucleon momentum p' = (\p+q\) > pp 

where pp n — [§7r 2 /On,p( r )] s , are the local Fermi mo- 
menta of neutrons and protons at the interaction point 
in the nucleus defined in terms of their respective nu- 
clear densities p n ,p(r). These constraints are incorpo- 
rated while performing the integration over the initial nu- 
cleon momentum in Eqn.(2) by replacing the energy con- 
serving 5- function in Eqn.(3) by — i-Zro£//v(g , o,<7) where 
Un(qo,q) is the Lindhard function corresponding to the 
particle hole(ph) excitations induced by the weak inter- 
action process through W exchange shown in Fig. 1(a). 
In the large mass limit of W boson i.e.Myy — * oo, this 
Fig. 1(a) is reduced to Fig. 1(b) for which the imaginary 
part of the Lindhard function i.e. ImUN(qo,q) is given 
by 
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FIG. 1: Diagrammatic representation of the neutrino self- 
energy diagram corresponding to the ph-excitation leading to 
v e + n — > e~ + p process in nuclei. In the large mass limit of 
the W boson (i.e.Mw —* oo) the diagram 1(a) is reduced to 
1(b) which is used to calculate \T\ 2 in Eqn.(5). 



1 M M 

ImU N (q ,q) = -- s£L<L [E Fl - A] with (8) 



2tt \q\ 



q 2 < 0, Ep 2 - q < E Fl and 



-qo + \q\\n- 



< E 



where Ep t and Ep 2 are the local Fermi energy of initial 
and final nuclcons and 



A = Max 



M n , Ep 2 - q Q , 



-ft + iau/i-^ 



The expression for the neutrino nuclear cross section 
a A (E v ), is then given by: 



<r A (E u ) = — 



E v Ei 



rPi /-l 
r 2 dr / pi dpi / d(cos9) 

J pi mi n J_X 

j2Y,\ T \ 2lmUN ^~ E ^ 3- ( 9 ) 



Moreover in the nucleus, the Q value of the nuclear reac- 
tion and the Coulomb distortion of the final lepton in the 
electromagnetic field of the final nucleus should be taken 
into account. This is done by modifying the energy con- 
serving (5-function 5{E V — Ei + E n — E p ) in Eqn.(3) to 
6(E„-Q - {Ei + V c (r)) + E n - E p ) where V c (r) is the 
Coulomb energy of the produced lepton in the field of 
final nucleus and is given by 



1 f 

VJr) = ZZ'a4ir(- / 
r Jo 



Z 



Jr % 

- . . (10) 

This amounts to evaluation of Lindhard function in 
Eqn.(8) at (q - (Q + V c (r)), q) instead of (q , q). The 
implementation of this modification requires a judicious 
choice of Q value for inclusive nuclear reactions in which 
many nuclear states are excited in iron. We have taken Q- 
value of 6.8MeV corresponding to the transition to lowest 
lying 1+ state in 56 Co for vi+ 56 Fe — > l~ + 56 Co* reaction 
and Q- value of 4.3MeV corresponding to the transition to 
lowest lying 1+ state in 56 Mn for vi+ b6 Fe -> l++ 56 Mn* 
reaction. 

The inclusion of V c (r) to modify energy and cor- 
responding momentum of the charged lepton in the 
Coulomb field of final nucleus in our model is equivalent 
to the treatment of Coulomb distortion effect in mod- 
ified effective momentum approximation(MEMA). This 
approximation has been used in other calculations of 
charged current neutrino reactions |36J and electron scat- 
tering at higher energies |37|. 

With these modifications, the final expression for the 
quasielastic inclusive production from iron nucleus is 
given by 



<J A {E V ) 



E v Ei 



rpi mam rl 

r 2 dr / p 2 dpi I d(cos9) 
^^|T| 2 /mlM<Zo-(Q + ^M), q](U) 



It is well known that weak transition strengths are mod- 
ified in the nuclear medium due to presence of strongly 
interacting nucleons. This modification of the weak tran- 
sitions strength in the nuclear medium is taken into ac- 
count by considering the propagation of particle- hole (ph) 
excitations in the medium. While propagating through 
the medium, the ph-excitations interact through the nu- 
cleon nucleon potential and create other particle hole and 
A-h excitations as shown in Fig. 2. The effect of these 
excitations are calculated in random phase approxima- 
tion which is described in Ref . [23l l32| . The effect of nu- 
clear medium on the renormalization of weak strengths is 
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(a) 



FIG. 2: Many body Feynman diagrams (drawn in the limit 
Mw — * oo) accounting for the medium polarization effects 
contributing to the process v £ + n — > e _ + p transitions. 



treated in a non-relativistic frame work. In leading order, 
the non-relativistic reduction of the weak hadronic cur- 
rent defined in Eqn.6, F2(q 2 ) term gives a spin-isospin 
transition operator t which is a transverse opera- 
tor while the Fa{c[ 2 ) term gives a spin-isospin transi- 
tion operator a ■ f which has a longitudinal as well as 
a transverse part. This representation of the transi- 
tion operators in the longitudinal and transverse parts 
is useful when summing the diagrams in Fig. (2) in ran- 
dom phase approximation(RPA) to calculate \T\ 2 . While 
the charge coupling remains unchanged due to nuclear 
medium effects, the terms proportional to i*f are affected 
by the transverse part of the nucleon nucleon potential 
while the terms proportional to F\ are affected by trans- 
verse as well as longitudinal parts. The effect is to re- 
place the terms like F 2 , F\, F 2 Fa etc. in the following 
manner 23. 321 



(F 2 ,F 2 F A )^(F 2 ,F 2 F A )- 



1 



U N V t \ 2 



3|i-c^v,r z\i-u N v t \ 



(12) 



where Vi and Vt are the longitudinal and transverse parts 
of the nucleon nucleon potential calculated with tt and p 
exchanges and modulated by Landau- Migdal parameter 
g' to take into account the short range correlation effects 
and are given by 



Vi{q) 



'A 2 



Vt(q) = hi 



in- 



Al-q 2 



Ap - rn 2 p 
A 2 - a 2 



with A T 



l.dGeV, C p = 2.0, A, 



2.5GeV, 



(13) 



and 



m p are the 



ion and rho meson masses and g' is taken 
to be 0.7|31|. The effect of Ah excitations are taken 
into account by including the Lindhard function Ua for 



the Ah excitations and replacing Upf by Un + Ua in 
Eqn.(12). The complete expression for Um and Ua used 
in our calculations are taken from 38] . The different 
couplings for N and A to the nucleon are incorporated 
in Um and Ua and then the same interaction strengths 
Vi and Vt are used for ph and Ah excit at ions |3fll] - |4*flj . 



B. Results 

We present the numerical results for the total cross 
section for the quasielastic processes vi{i>i) + 56 Fe — > 
+ 56 Co*( 5& Mn k ) as a function of energy for neu- 
trino and anti neutrino reactions on iron in the energy 
region relevant to the fully contained events of atmo- 
spheric neutrinos i.e. E v < 3GeV. The cross sections 
have been calculated using Eqn.(ll) with the nuclear- 
density p(r) given by a two parameter Fermi density |4lj: 



with c = 3.971/m, z = 0.5935/m, 



{A 



jp-pir) and p p (r) = fp(r) 



p{r) = 
Pn{r) = 

In Fig. 3 we show the numerical results of u(E) vs E, 
for all flavors of neutrinos i.e. f M , v^, v e and v e . The 
reduction due to nuclear effects is large at lower ener- 
gies but becomes small at higher energies. The energy 
dependence of the cross sections for muon and electron 
type neutrinos are similar except for the threshold effects 
which are seen only at low energies (E v < 50QMeV). 
This reduction in a is due to Pauli blocking as well as due 
to the weak renormalization of transition strengths which 
have been separately shown in Fig. 4(a) and Fig. 4(b) 
for neutrinos and antineutrinos, where we also show 
the results in the Fermi gas model given by Llewellyn 
Smith[2(j. We plot in Fig. 4(a) and Fig. 4(b) for neutrinos 



and antineutrinos. the reduction factor R 



ar{E) 



vs E where cr n udear{E) is the cross section per neu- 
tron(proton) for neutrino (antineutrino) reactions in the 
nuclear medium. The solid lines show the reduction fac- 
tor R when only the Pauli suppression is taken into ac- 
count through the imaginary part of the Lindhard func- 
tion given in Eqn.8. This is similar to the results of 
Llewellyn Smith 20] in Fermi gas model shown by dashed 
lines. In this model the total cross section a is calculated 
by using the formula a = J dq 2 R(q 2 ){j^i) free where 
R(q 2 ) describes the reduction in the cross section calcu- 
lated in the Fermi gas model and includes the effect of 
Pauli suppression onlv|20|. However, in our model we get 
further reduction due to renormalization of weak transi- 
tion strengths in the nuclear medium when the effects 
of Fig. 2(a) and Fig. 2(b) are included. These are shown 
by dotted lines in Figs. 4(a) and 4(b). We find that the 
reduction at higher energies (E > lGeV) is around 20% 
for neutrinos and 40% for antineutrinos. It is worth not- 
ing that the energy dependence of the reduction due to 
nuclear medium effects is different for neutrinos and an- 
tineutrinos. This is due to the different renormalization 
of various terms like F\, F 2 Fa and F 2 in |T| 2 which 
enter in different combinations for neutrino and antincu- 
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trino reactions. The results for v e and v e cross section 
are respectively similar to and reactions except for 
the threshold effects and are not shown here. 

In Fig. 5 and Fig. 6 we compare our results for a{E) with 
the results of some earlier experiments which contain 
nuclear targ ets like Carbo nlip . Freonji3|-0|, Freon- 
Propane|45j and Aluminum |46j|. where the experimental 
results for the deuteron targets H3 are not included as 
they are not subject to the various nuclear effects dis- 
cussed here. It should be kept in mind that the nu- 
clear targets considered here (except for Br in Freon) 
are lighter than Fe. Therefore, the reduction in the total 
cross section due to nuclear effects will be slightly over- 
estimated. For example, for energies E v > lGeV the re- 
duction in neutrino(antineutrino) cross section in case of 
56 Fe is 5% more than the reduction in case of 12 C0 . In 
comparision to the neutrino (antincutrino) nuclear cross 
sections as obtained in the Fermi gas model of Llewellyn 
Smith 20] (shown by dashed lines in Fig. 5 and Fig. 6) we 
get a smaller result for these cross sections. This reduc- 
tion in the total cross section leads to an improved agree- 
ment with the experimental results as compared to the 
Fermi gas model results specially for antineutrino reac- 
tions(Fig.6). It should be emphasized that the Fermi gas 
model has no specific mechanism to estimate the renor- 
malization of weak transition strengths in nuclei while in 
our model this is incorporated by taking into account the 
RPA correlations. 

In Fig. 7(a) and Fig. 7(b), we show the nuclear medium 
effects on the momentum and angular distributions i.e. 
^ and d fj s6[ of leptons produced in and reactions. 
We find a large suppression in the results specially in 
the peak region of momentum and angular distributions. 
Quantitatively similar results are obtained for the case of 
v e and v e reactions and are not shown here. 

We have also studied the effect of Coulomb distortion 
in the momentum distribution of leptons but find no sub- 
stantial effect around E v — l.OGeV. These are found to 
affect the results only at low energies i.e. E v < 500MeV 
where the peak is slightly shifted to lower momentum as 
shown in Figs. 8(a) and 8(b). 



III. INELASTIC PRODUCTION OF LEPTONS 

The inelastic production process of leptons is the pro- 
cess in which the production of leptons is accompanied by 
one pion (or more pions). There are many calculations of 
one pion production by neutrinos from free nucleons|49j 
but there are only a few calculations which discuss the 
nuclear effects in these processes l50|- [521. In this sec- 
tion we follow the method of Ref.[&I| to estimate the 
nuclear effects and nuclear model dependence of inelastic 
production cross section of leptons induced by neutrinos 
from iron nuclei. The calculations are done assuming 
A-dominance of one pion production because the contri- 
bution of higher resonances in the energy region of at- 
mospheric neutrinos leading to fully contained events is 




400 600 
E v (MeV) 



FIG. 3: Total quasielastic cross sections in the present model 
for the neutrino and antineutrino reactions in 56 Fe are shown 
by solid line for v e , dashed line for v e , dotted line for and 
dashed-dotted line for reactions. 



expected to be small. 



A. Formalism 

The matrix element for neutrino production reaction 
of A on proton targets leading to one pion events i.e. 



v l (k)+p( P )^r(k , ) + A++(p') 



(14) 



is given by Eqn.5 where J M now defines the matrix ele- 
ment of the transition hadronic current between N and 
A states. The most general form of Jq C is written as[5lj : 



J£ = Mp') 



g 3 v (g 2 ) 

M 



+ 



^p-{g a ^ ■ p - ?V) + ^§%V j 



M 



M 2 



m yy H q ' ' M 2 



+ C£{q 2 )g^ 



M 2 



-q <T 



*(P) 



(15) 



where ip a (p') and u(p) are the Rarita Schwinger and 
Dirac spinors for A and nucleon of momenta p' and p re- 
spectively, q(= p' —p = k — k') is the momentum transfer 
and CY (i=3-6) are vector and C/ l (i=3-6) are axial vector 
transition form factors. The vector form factors C^(i=3- 
6) are determined by using the conserved vector cur- 
rent(CVC) hypothesis which gives Cg / (g 2 )=0 and relates 
Cj (i=3,4,5) to the electromagnetic form factors which 
are determined from photoproduction and electroproduc- 
tion of A's. Using the analysis of these experiments!^- 
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FIG. 4: Ratio of the total cross section to the free neutrino nucleon cross section for the reactions (a) + n —* + p (b) 
+ P — » M + + n in iron nuclei in the present model with Pauli suppression (solid line) and with nuclear effects (dotted line) 
and in the Fermi gas model (dashed line) |20| . 




E v (GeV) 



FIG. 5: Neutrino quasielastic total cross section per nucleon in iron for + n — > p + n reaction. The data are from 
LSND|!2 (Ellipse), Bonnetti et al. HI (squares), SKAT collab.0 (Triangle Down), Pohl et alj|l (Circle) and Belikov et 
al. 46] (Triangle Up). The dashed line is the result of the cross section in the Fermi gas model 20] and solid line is the re- 
sult using the present model with nuclear effects. 



0.8 




12 3 4 

E v (GeV) 



FIG. 6: Anti-neutrino quasielastic total cross section per nucleon in iron for + p — > n + /x + reaction. The data are from 
Bonnetti et al-Hl, SKAT collabn.|3 (Triangle Down), Pohl et al.|H] (Circle), Belikov et al. Hi (Triangle Up). The dashed line 
is the result of the cross section in the Fermi gas model|2(J and solid line is the result using the present model with nuclear 
effects. 
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FIG. 7: The results for the differential cross section Fig. (a) j£- vs p; and Fig.(b) d ^ s6 vs cosft at E = l.OGel/ in the present 
model with Pauli suppression (solid line for and dashed-dotted line for z/ M ) and with nuclear effects (dashed line for and 
dotted line for Vp). 




FIG. 8: Differential scattering cross section vspi for E v = 500MeV, with Coulomb effect (solid line for v% and dashed-dotted 
line for vi) and without Coulomb effect (dashed line for vi and dotted line for v{) Coulomb effect. Fig. (a) is for electron type 
and Fig.(b) for muon type neutrinos. 



[53j we take for the vector form factors 

M 



cY = o, cr = 



v 



M A 



-C 



v 

3 ' 



and 



CSV) = 



2.05 



[ Ml I 



Ml = OMGeV 2 



(16) 



The axial vector form factor Cg(q 2 ) is related to C$ A (q 2 ) 
using PC AC and is given by 



cW) = ckV)- m2 o 



M 2 



The remaining axial vector form factor C^ 3 4 5 (g 2 ) 
are taken from the experimental analysis of the neu- 
trino experiments producing A's in proton and deuteron 
targets [54| -|55|. These form factors are not uniquely de- 
termined but the following parameterizations give a sat- 



isfactory fit to the data. 



Ct 3A , 5 (q 2 ) = Cf(0) 



Ma 



(18) 

with Cf(Q) = 0,(7^(0) = -0.3,C^(0) = 1.2, a 4 = a 5 = 
-1.21, b 4 = b 5 = 2GeV 2 , M A = l.28GeV. Using the 
hadronic current given in Eqn.(15), the energy spectrum 
of the outgoing leptons is given by 



d 2 a 



1 



1 k! 



2 



(17) dE k -dn k , Sir 3 MM' E„ ( W - M') 2 



r 2 (w) v 

4. 



L P lL/ 



(19) 



where W = J (p + q) 2 and M' is mass of A, 

Lfi V = k^k v + k^k u 9fii/k ■ k + ic^yapk k ^ , 



8 



and is calculated with the use of spin | projection oper- 
ator defined as 



and given by: 



2M> \ 9 3 M' 2 
Ip"*7"-P"'7p 1 



M' 



(20) 



In Eqn.(19), the decay width T is taken to be an energy 
dependent P-wave decay width given by 

T (W) = ±- ( ^ cm \ 3 Q(W - M - m.) (21) 

D7T \ m,r / W 

where 

VQE 2 - m 2 - Wf - 4m 2 M 2 

and M is the mass of nucleon. The step function 9 de- 
notes the fact that the width is zero for the invariant 
masses below the Nir threshold. |q cm | is the pion mo- 
mentum in the rest frame of the resonance. When the re- 
action 14) takes place in the nucleus, the neutrino inter- 
acts with the nucleon moving inside the nucleus of density 
p(r) with its corresponding momentum p constrained to 
be below its Fermi momentum. The produced A's have 
no such constraints on their momentum. These A's de- 
cay through various decay channels in the medium. The 
most prominent decay mode is A — > Nit which produces 
pions. This decay mode in the nuclear medium is slightly 
inhibited due to Pauli blocking of the final nucleon mo- 
mentum modifying the decay width T used in Eqn.(21). 
This modification of T due to Pauli blocking of nucleus 
has been studied in detail in electromagnetic and strong 
interactions 56J. The modified A decay width i.e. T is 
written as 56J: 



1 {Una\ M\ 

Hem I 



w 



6lT \ TOtt 



r = — j j '--/•.:/, 7 . /-.a./.-aiBMI .\/ m n ) 

(22) 

where F(kp, E A , &a) is the Pauli correction factor given 
by 

kAKm\ + E A E' — EpW 
F(k F ,E A: k A )= AlHcml - , , Pc 7 — (23) 



2fc A |q' c 



where fcj? is the Fermi momentum, Ef = y M 2 



is the A momentum and E A = WW 

Moreover, in the nuclear medium there are addi- 
tional decay channels now open due to two body and 
three body absorption processes like AiV — * NN and 
ANN — > NNN through which A's disappear in the nu- 
clear medium without producing a pion while a two body 



A absorption process like AN — ► ttNN gives rise to some 
more pions. These nuclear medium effects on A propa- 
gation are included by using a A propagator in which 
the A propagator is written in terms of A self energy 
Sa- This is done by using a modified mass and width 
of A in nuclear medium i.e. M A — > M\ + ReT, A and 
f — > f — JtoEa. There are many calculations of A self 
energy S a in the nuclear medium |56|-|59| and we use the 
results of |56| , where the density dependence of real and 
imaginary parts of Sa are parametrized in the following 
form: 

i?e£ A = 40— MeV and 
Po 



— JtoEa = Cq 



,P / \ 7 

P \ , ^ IP 



PO 



C az\jJ (24) 



In Eqn.24, the term with Cq accounts for the AN — > 
ttNN process, the term with Cai for two-body absorp- 
tion process AiV — > NN and the term with Ca3 for 
three-body absorption process ANN — > NNN. The 
coefficients Cq, Cai, Ca3, ot, (3 and 7 (7 = 2/3) are 
parametrized in the range 80 < < 320MeV (where 
T v is the pion kinetic energy) as |56J 

Ci(T n ) = ax 2 +bx + c, for x = — (25) 

The values of coefficients a, b and c are given in Table-I. 
taken from ref. 561. 



TABLE I: Coefficients of Eqn.(25) for an analytical interpo- 
lation of ImTfA 





C Q (MeV) 


CA 2 (MeV) 


C A3 (MeV) 


a 


P 




-5.19 


1.06 


-13.46 


0.382 


-0.038 


b 


15.35 


-6.64 


46.17 


-1.322 


0.204 


c 


2.06 


22.66 


-20.34 


1.466 


0.613 



With these modifications, which incorporate the vari- 
ous nuclear medium effects on A propagation, the cross 
section is now written as 



dv dk' 1 

8tt 3 E v Ei MM' 



(W 



M' - i?e£ A ) 2 + (57 - Itu^a) 2 
1 



:Pn(r) 



(26) 



The factor | in front of p n comes due to suppression of 
charged pion production from neutron targets i.e. v\ + 
n — > l~ + A + — > l~~ + n + n + , as compared to the charged 
pion production from the proton target i.e.i^ +p — ► l~ + 
A ++ — > l~ +p+n + , in the nucleus. In case of antincutrino 
reactions p p + \p n is replaced by p n + \p p - 
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B. Results 

In this section we present results of inelastic lepton 
production cross section due to Ah excitations in iron 
induced by the charged current neutrino interactions us- 
ing Eqn.(26). In Fig. 9(a) and Fig. 9(b), we show the 
results <j{E v ) ~ E v for v e (v e ) and v^iy^ for the in- 
elastic production of lepton accompanied with A and 
compare this with the cross section for the quasielas- 
tic production of leptons discussed in section-II. We see 
that for E v « lAGeV the lepton production cross sec- 
tion through quasielastic and inelastic production pro- 
cesses are comparable. For energies E v < lAGeV where 
the atmospheric neutrino energies are important for fully 
contained events the major contribution comes from the 
quasielastic events. 

The effect of nuclear effects on the A production are 
shown in Fig. 10 for and Dp reactions. The results for 
the v e and v e are similar to Fig. 10 and are not shown 
here. We see that the nuclear medium effects reduce the 
A production cross section by 5-10%. 

In Fig. 11 (a) and Fig. 11(b), we show the momentum 
distribution da / dp and angular distribution do/dcosO for 
muon type neutrinos where we also show the effects of 
nuclear effects. The nuclear medium effects reduce the 
cross sections mainly in the peak region of the momen- 
tum and angular distributions for muons. The results for 
the momentum distribution and angular distribution for 
electrons is similar to the muon distributions. We show 
our final results on momentum and angular distributions 
for all charged leptons i.e. e - , e + , fjT and fi + produced in 
neutrino and antineutrino reactions with u e , v^, v e and 
Up for E v = lGeV in Fig.l2(a) and Fig.l2(b) in 56 Fe 
nuclei. 

Here we will like to make some comments about the 
lepton events produced through the A excitation. The 
A excitation process gives rise to leptons accompanied by 
one pion events produced by A — > Ntt and AiV — > irNN 
processes in the nuclear medium. The pions produced 
through these processes will undergo secondary nuclear 
interactions like multiple scattering and possible absorp- 
tion in iron nuclei while passing through the nucleus and 
an appropriate model has to be used for their description. 
Models developed by Salcedo et al.|!3 and Paschos et 
have studied these effects but we do not consider 
these effects here as they are beyond the scope of this 
paper. The A excitation process in the nuclear medium 
also gives rise to quasielastic like events through two body 
and three body absorption processes like A7V — > NN and 
AiV — >■ ANN where only leptons are present in the final 
states. The quasielastic-like lepton events are discussed 
by Kim et al.[21j but no quantitative estimates have been 
made. We have in an earlier paper [SlJ discussed this 
process only qualitatively but make a quantitative esti- 
mate of these events in this paper. We find that around 
E v = lGeV the contribution of these quasielastic like 
events is about 10 — 12% but its effect on the flux aver- 
aged production of leptons for atmospheric neutrinos is 



not very significant. This is discussed in some detail in 
the next section. 



IV. LEPTON PRODUCTION BY 
ATMOSPHERIC NEUTRINOS 

The energy dependences of the quasielastic and inelas- 
tic lepton production cross sections described in sections 
II and III have been used to calculate the lepton produc- 
tion by atmospheric neutrinos after averaging over the 
neutrino flux corresponding to the two sites of Soudan 
and Gransasso, where iron based detectors are being 
used. There are quite a few calculations of atmospheric 
neutrino fluxes at these two sites. We use the angle av- 
eraged fluxes calculated by Honda et al.Q and Barr 
et al-ITH for the Soudan site and the fluxes of Barr et 
al.|l5j and Plyaskin 16J for the Gransasso site to calcu- 
late the flux averaged cross section < a > and also the 
momentum and the angular distributions < 4 s - > and 

& dpi 

< rfc pgg > for leptons produced by v e , u e , and v^. 



A. Flux averaged Momentum and Angular 
distributions 

In this section we present the numerical results for the 
flux averaged momentum distributions < ^ > and an- 
gular distributions < d f° s8 > for various leptons pro- 
duced from v e , 24, i/p, at the atmospheric neutrino 
sites of Soudan and Gransasso. These leptons are pro- 
duced through the quasielastic as well as inelastic pro- 
cesses. The quasielastic production has been discussed 
in section-II and the inelastic production has been dis- 
cussed in section-Ill. The inelastic production of lep- 
tons is accompanied by pions. However, in the nu- 
clear medium where the production of A is followed by 
AN -> NN, ANN -> NNN absorption processes, the 
leptons are produced without pions. These are quasielas- 
tic like events. We show the momentum distribution of 
the leptons produced in iron nuclei corresponding to one 
pion and quasielastic like events for atmospheric neutri- 
nos relevant to Soudan and Gransasso sites. We show 
it for and for the Soudan site in Figs. 13(a) and 
13(b) and for the Gransasso site in Figs. 14(a) and 14(b) 
corresponding to the flux of Barr et al.jl^. We see that 
at both sites, the production cross section of quasielastic- 
like events is quite small. 

We now present our final results for the momentum 
distribution and angular distribution for various lep- 
tons produced by atmospheric neutrinos at Soudan and 
Gransasso sites in Figs. 15-18. In these figures sep- 
arate contributions from the quasielastic and inelastic 
processes to the momentum and angular distributions 
of charged leptons are shown explicitly. The quasielas- 
tic events are those where only a charge lepton is pro- 
duced in the final state either by a quasielastic process 
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described in section-II or by an inelastic process of A- 
production followed by its subsequent absorption in the 
nuclear medium as described in section-Ill. The inelastic 
events are those events in which a charged lepton in the 
final state is accompanied by a charged pion as a decay 
product of deltas excited in the nuclear medium. 

We show the results for the momentum distribution 
< da /dpi >~ pi(l — e~, e + , for the atmospheric 

neutrino fluxes of Barr et al.[15| at Soudan site in Fig. 15 
and at Gransasso site in Fig. 16. We see that in all cases 
the major contributions to the charged lepton production 
comes from the quasielastic processes induced by neu- 
trinos (solid line). The contribution from the quasielas- 
tic processes induced by antineutrinos (dotted line) and 
inelastic processes induced by neutrinos (dashed line) is 
small and is around 20% in the peak region. The contri- 
bution due to inelastic processes induced by antineutrinos 
is very small over the whole region (dashed-dotted lines). 

The peak in quasielastic reactions occurs around 
pi w 200MeV. The peaks in the inelastic and 
quasielastic v„ reactions are slightly shifted towards lower 
energies. The momentum distribution of the leptons for 
the quasielastic reaction is peaked more sharply than the 
momentum distribution of the inelastic reaction. 

We have also presented the numerical results for an- 
gular distributions of leptons < , da „ > vs cosO for 

° * dcos&i 

the atmospheric neutrino fluxes of Barr et al.[15| for the 
Soudan site in Fig. 17 and for the Gransasso site in Fig. 18. 
The inelastic lepton production accompanied by pions 
(dashed line for /i~(e~) production and dashed-dotted 
line for /i + (e + ) production) and the quasielastic lepton 
production events (solid line for /i~(e~) production and 
dotted line for /i + (e + ) production) have been explicitly 
shown in these figures. We see from these figures that the 
lepton production cross sections are forward peaked in 
all cases. The inelastic distribution due to pion produc- 
tion are slightly more forward peaked than the quasielas- 
tic distribution. The contribution of the inelastic lepton 
events are small compared to the quasielastic events and 
the angular distribution of the flux averaged cross sec- 
tions are dominated by the quasielastic events. 

At Soudan site, we have also studied the momentum 
and angular distributions for the atmospheric neutrino 
fluxes of Honda et al. 01 • We find that the momentum 
and angular distributions for the production of muons are 
similar to the distributions obtained for the flux of Barr 
et al.[15|. In the case of electron production, the use of 
the flux of Honda et al. gives a slightly smaller value 
for < da/dpi > and < d £° sBi > for electrons as compared 
to the flux of Barr et al.flfij. 

Similarly at Gransasso site, the momentum and an- 
gular distributions for the atmospheric neutrino fluxes 
of Plyaskin^^ have also been studied. Here, also we 
find that the momentum and angular distributions for 
the production of muons are similar to the distributions 
obtained for the flux of Barr et al. |l5j . In the case of elec- 
tron production, the use of the flux of Plyaskin|l6|, gives 
a slightly smaller value for < da /dpi > and < d f° sg > 



for electrons as compared to the flux of Barr et al.ll5|. 

We further find that for the flux of Barr et al.jlq. 
the lepton production cross section at the Soudan site is 
slightly larger than the Gransasso site for muons as well 
as for electrons. 



B. Flux averaged Total cross sections and lepton 
yields 

The total cross sections for the production of leptons 
and its energy dependence have been discussed in section 
II and section III for quasielastic and inelastic reactions. 
In this section we calculate the lepton yields Yi for lepton 
of flavor I which we define as 

Yi = J <?{E Vl ) dE Vl 

where, is the atmospheric neutrino flux of vi and 
a(E Vl ) is the total cross section for neutrino vi of energy 
E Vl . We calculate this yield separately for the quasielas- 
tic and inelastic events. We define a relative yield of 
muon over electron type events by R as R = R^/ e = 

Y + y- f° r quasielastic and inelastic events and present 
our results in Table-II. We study the nuclear model de- 
pendence as well as the flux dependence of the relative 
yield R. 

The results for R are presented separately for 
quasielastic events, inelastic events and the total events 
in Tablc-II. For quasielastic events vi(yi) + 56 Fe — > 
+ X, the results are presented for the case of free 
nucleon by Rfn, for the nuclear case with Fermi gas 
model description of Llewellyn Smith |2(J by Rfg and 
for the case of nuclear effects within our model by Rnm- 
We see that there is practically no nuclear model depen- 
dence on the value of R (compare the values of Rnm, 
Rfg and Rfn for the same fluxes at each site). This 
is also true for the inelastic production of leptons i.e. 
^(^)+ 56 Fe^ r(l+)+w + (tt") + X for which the ratio 
for free nucleon case (denoted by Raf) and the ratio for 
the nuclear case in our model (denoted by Ran) are pre- 
sented in Table II. It is, therefore, concluded that there is 
no appreciable nuclear model dependence on the ratio of 
total lepton yields for the production of muons and elec- 
trons (compare the values of Rf and Rn, where Rf shows 
the ratio of total lepton yields for muon and electron for 
the case of free nucleon and Rn shows the ratio of total 
yield for muon and electron for the case of nucleon in the 
nuclear medium). 

However, there is some dependence of the ratio R on 
the atmospheric neutrino fluxes. The flux dependence of 
R can be readily seen from Table-II, for the two sites of 
Soudan and Gransasso. At the Gransasso site, we see 
that there is 4-5% difference in the value of Rn for the 
total lepton yields for the fluxes of Barr et al.[l5| and 
Plvaskin[l6j. At Soudan site, the results for the fluxes of 
Honda et al.^4| and Barr et al.E^] are within 4-5% but 
the flux calculation of Plyaskin 16j gives a result which 
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is about 10-11% smaller than the results of Honda et 
al.[l4| and 7-8% smaller than the results of Barr et al.[l5j. 
The flux dependence is mainly due to the quasielastic 
events. This should be kept in mind while using the flux 
of Plyaskin[3 for making any analysis of the neutrino 
oscillation experiments. 

In Table-Ill, we present a quantitative estimate of the 
relative yield of inelastic events 77 defined by r; = -^7-, 

where Y/ A = Y/ A + Yj A is the lepton yield due to the 
inelastic events and Yi is the total lepton yield due to the 
quasielastic and inelastic events i.e.Y/ = Yi q e ' + Yj q e ' + 
Yi A + Yj A . The relative yield for the case of free nucleon 
is shown by ri(F) and for the case with the nuclear ef- 
fects in our model is shown by 77(A). We see that for 
free nucleons, the relative yield of the inelastic events 
due to A excitation is in the range of 12-15% for vari- 
ous fluxes at the two sites. The ratio is approximately 
same for electrons and muons. When the nuclear effects 
are taken into account this becomes 19-22%. This is due 
to different nature of the effect of nuclear structure on 
the quasielastic and inelastic production cross scetions 
which gives a larger reduction in the cross section for the 
quasielastic case as compared to the inelastic case. This 
quantitative estimate of 77 (A) in iron may be compared 
with the results in oxygen where the experimental results 
at Kamiokande give a value of 18%|60|. 

V. CONCLUSIONS 

We have studied the charged lepton production in iron 
induced by atmospheric neutrinos at the experimental 
sites of Soudan and Gransasso. The energy dependence 
of the total cross sections for the quasielastic and inelastic 
processes have been calculated in a nuclear model which 
takes into account the effect of Pauli principle, Fermi 
motion effects and the renormalization of weak transi- 
tion strengths in nuclear medium. The inelastic process 
has been studied through the A dominance model which 
incorporates the modification of mass and width of A res- 
onance in the nuclear medium. The numerical results for 
the momentum and angular distributions of the charged 
lepton production cross section have been presented for 
muons and electrons. The relative yield of muon to elec- 
tron production has been studied. In addition to the nu- 
clear model dependence, the flux dependence of the total 
yields, momentum distribution 4 s - and , da „ have been 

J ' dpi acos&i 

also studied. In the following we conclude this paper by 
summarizing our main results: 

1. There is a large reduction due to nuclear effects in the 



total cross section for quasielastic production cross sec- 
tion specially at lower energies(40-50% around 200MeV) 
and the reduction becomes smaller at higher energies (15- 
20% around 500MeV ). 

2. For quasielastic reactions we find a larger reduction 
in the total cross section as compared to the Fermi gas 
model. The energy dependence of this reduction in cross 
section at low energies (E < 500MeV) is found to be 
different for neutrino and antineutrino reactions. 

3. The inelastic production of leptons where a charged 
lepton is accompanied by a pion becomes comparable to 
the quasielastic production of leptons for E v ~ lAGeV 
and increases with increase of energy. The effect of nu- 
clear medium on the inelastic production cross section 
(in absence of pion re-scattering effects) in not too large 
(~ 10%). 

4. For quasielastic reactions the effect of Coulomb dis- 
tortion of the final lepton in the total cross section is 
small except at very low energies (E < 500MeV) and 
becomes negligible when averaged over the flux of atmo- 
spheric neutrinos. 

5. The flux averaged momentum distribution of leptons 
produced by atmospheric neutrinos is peaked around the 
momentum pi ~ 200MeV for electrons and muons. The 
peak for the quasielastic production is sharper than the 
inelastic production. The inelastic production of leptons 
contributes about 20% to the total production of leptons. 

6. The flux averaged angular distribution of leptons for 
atmospheric neutrinos for quasielastic as well as inelastic 
production is sharply peaked in the forward direction. 

7. There is a very little flux dependence on the relative 
yield of muons and electrons at the site of Gransasso. 
However, at the Soudan site, the atmospheric flux as de- 
termined by Plyaskin^ij gives a value of relative yield 
which is smaller than the relative yield using the flux of 
Honda et al.^3| an d Barr et al.|l5j. 

8. The nuclear model dependence of the relative yield of 
muons to electrons is negligible, even though the individ- 
ual yields for muons and electrons are reduced with the 
inclusion of nuclear effects specially for the quasielastic 
production. 
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vi + N — » A + l + reactions (shown by dashed-dotted line) in 56 Fe and compared with the corresponding quasielastic cross 
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Fig. (a) is for electron type and Fig.(b) is for muon type neutrinos. 
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FIG. 11: Differential scattering cross section Fig. (a) vs pi and Fig.(b) d fj T s g [ vs cosOi at E v — l.OGeV for the reaction 
viipi) + N — > A + in 56 Fe with (dashed line for dotted for P M ) and without nuclear effects(solid line is the result for 

reaction and dashed-dotted for P^). 
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FIG. 12: The results for the differential cross section Fig. (a) -jj^- vs pi and Fig.(b) d ff s0l vs cosOi with nuclear effects at 
E — l.OGeV for the reaction vi(p{) + N — > A + in 56 Fe. Solid line is the result for u e , dashed for v^, dotted for i> e and 

dashed-dotted for P M reactions. 




FIG. 13: Flux averaged differential cross section < > vs pi at the Soudan site simulated by Barr et al.[T^|. Fig. (a) for 

+ N — > A + and Fig.(b) for i/ M + iV — > A + M + reactions in 56 Fe. The solid line is for total A events, dashed line is for 
one 7T production events and dotted line is for quasielastic like events(due to A absorption in the nuclear medium) 
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FIG. 14: Flux averaged differential cross section < > vs pi at the Gransasso site simulated by Barr et al.|l5|. Fig. (a) for 

U/j, + N — > A + /i~ and Fig.(b) for P M + N — ► A + M + reaction in 56 .Fe. The solid line is for total A events, dashed line is for 
one 7r production events and dotted line is for quasielastic like events(due to A absorption in the nuclear medium) 
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FIG. 15: Flux averaged differential cross section < > vs pi at the Soudan site simulated by Barr et al.[l5T|. The solid line is 

for total quasielastic events and dashed line is the one 7r production events for z/ M reaction in 56 Fe. The corresponding results 
for Tip are shown by dotted and dashed-dotted lines, Fig. (a) for muon type and Fig.(b) for electron type neutrinos. 




FIG. 16: Flux averaged differential cross section < > vs pi at the Gransasso site simulated by Barr et al.[l3|. The solid 

line is for total quasielastic events and dashed line is the one 7r production events for reaction in 56 Fe. The corresponding 
results for are shown by dotted and dashed-dotted lines, Fig. (a) for muon type and Fig.(b) for electron type neutrinos. 
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FIG. 17: Flux averaged differential cross section < d ^° aBl > vs cosOi at the Soudan site simulated by Barr et al.jla]. The solid 
line is for total quasielastic events and dashed line is the one n production events for reaction in 56 Fe. The corresponding 
results for are shown by dotted and dashed-dotted lines, Fig. (a) for muon type and Fig.(b) for electron type neutrinos. 
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FIG. 18: Flux averaged differential cross section < de ^ 9 > vs cosOi at the Gransasso site simulated by Barr et al.|l5l|. The solid 

line is for total quasielastic events and dashed line is the one n production events for reaction in 56 Fe. The corresponding 
results for are shown by dotted and dashed-dotted lines, Fig. (a) for muon type and Fig.(b) for electron type neutrinos. 
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